Might we eventually understand the origin of the dark matter 

velocity anisotropy? 



oo 
O 
O 
(N 

O 

Q 

Oh 

6 
n 



> 
OO 

o 

00 

o 

> 

X 
S3 



Steen H. Hansen 

Dark Cosmology Centre, Niels Bohr Institute, University of Copenhagen, 
Juliane Maries Vej 30, 2100 Copenhagen, Denmark 
hansenOdark-cosmology . dk 



ABSTRACT 

The density profile of simulated dark matter structures is fairly well-established, and several 
explanations for its characteristics have been put forward. In contrast, the radial variation of 
the velocity anisotropy has still not been explained. We suggest a very simple origin, based on 
the shapes of the velocity distributions functions, which are shown to differ between the radial 
and tangential directions. This allows us to derive a radial variation of the anisotropy profile 
which is in good agreement with both simulations and observations. One of the consequences of 
this suggestion is that the velocity anisotropy is entirely determined once the density profile is 
known. We demonstrate how this explains the origin of the 7-/? relation, which is the connection 
between the slope of the density profile and the velocity anisotropy. These findings provide us 
with a powerful tool, which allows us to close the Jeans equations. 

Subject headings: 



1. Introduction 

The natural outcome of cosmological structure 
formation theory is equilibrated dark matter (DM) 
structures. According to numerical simulations, 
the mass density profile, p(r), of these structures 
changes from something with a fairly shallow pro- 
file in the central region, 7 = d\np/d\nr ~ —1 (or 
maybe zero) , to something steeper in the outer re- 
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Henriksen 2007 20081, varying from central cores 



to cusps. 

The second natural quantity to consider (af- 
ter the density profile) is the velocity anisotropy, 
which is defined through 



structures, like galaxy clusters, there appears to 
be fair agreement between the numerical pre- 
dictions and observations concerning the central 
steepness ( Pointecouteau et al.||2005 Sand et al. 
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where of 



and of are the 1-dimensional tangen- 



Vikhlinin et al. 20061, however, for smaller struc- 



tures, like galaxies or dwarf galaxies, observations 



tend to indicate central cores (Salucci et al. 2003 



tial and radial velocity dispersions (Binney & 
Tremaine 1987). If most dark matter particles in 
an equilibrated structure were purely on radial or- 
bits, then (3 could be as large as 1, and for mainly 
tangential orbits j3 could be arbitrarily large and 
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negative. Since dark matter is collision-less (3 does 
not have to be zero, and it could in principle even 
vary as a function of radius. 

Numerical N-body simulations of collision-less 
dark matter particles show that the dark matter 
velocity anisotropy is indeed radially varying, and 
that (3 goes from roughly zero in the central region, 



to 0.5 towards the outer region (Cole & Lacey 



1996 Carlberg et al. 1997). Only very recently has 



this velocity anisotropy been measured to be non- 
zero in galaxy clusters ( Hansen fc Piffaretti|2007 1 , 
and it has even been observed to be increasing 



as a function of radius (Host et al. 20091, in ex- 



cellent agreement with the numerical predictions. 
For smaller structures, like our own galaxy, this 
has not been observed yet. In principle (3 of our 
Galaxy can be measured in an underground direc- 
tional sensitive detector, however, it will require a 



large dedicated experimental programme (Host & 
Hansen||2007 1 . Very little theoretical understand- 
ing of the origin of this velocity anisotropy ex- 
ists, and to my knowledge no successful derivation 



of it has been published (see, however, Hansen & 
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2008)). We will in this paper present an 



attempt towards deriving (3. 

2. Decomposition 

When analyzing the outcome of a numerically 
simulated dark matter structure one traditionally 
divides the equilibrated structure in bins (shells) 
in radius, or in potential energy. For spherical 
structures there is naturally no difference. We can 
now consider all the particles in a given radial bin, 
and calculate properties like average density, an- 
gular momentum, velocity anisotropy etc. In order 
to do this, we must decompose the velocity of each 
particle into the radial component, and the two 
tangential components. The two tangential com- 
ponents can for instance be separated according 
to the total angular momentum of all the particles 
in the bin. 

By summing over all the particles in the radial 
(or potential) bin, we thus get the velocity distri- 
bution function (VDF), which for a gas would have 
been a Gaussian represented by the local gas tem- 
perature, f(v) ~ exp(—E/T). We are here dis- 
cussing the 1-dimensional VDF (i.e. the one where 
the two other velocities are integrated over), and 



we are not assuming that the radial and tangen- 
tial VDF's are independent. Naturally, since dark 
matter particles are not collisional, the concept of 
temperature is not well defined for them. In nu- 
merically simulated structures one observes that 
the radial VDF is symmetric (with respect to par- 
ticles moving in or out of the structure), and also 
the non-rotational part of the tangential VDF is 
symmetric. The asymmetry of the rotational part 
of the tangential VDF was discussed in Schmidt 
et all (120081) for the DM particles. For the struc- 



tures with very little rotation, the two tangential 
VDF's are virtually identical. To avoid any com- 
plications from the total angular momentum we 
will hereafter only discuss the two symmetric 1- 
dimcnsional VDF's. For any given radial bin in 
a given DM structure, the shape of the VDF only 
depends on the moment aneous distribution of par- 
ticles (which should be virtually time-independent 
for equilibrated structures) , and is independent of 
the method by which the structure is selected. 

When analyzing dark matter structures result- 
ing from cosmological simulations, we find that 
the shape of the radial VDF changes as a func- 
tion of radius (Wojtak et al. 2005| Hansen et al 



2006] IFaltenbacher fe Diemand||2006| |Fairbairn~fe" 
Schwetz|2008 ). In particular, the bins in the inner 



region tend to have long tails (more particles at 
high velocity compared to a Gaussian), whereas 
bins at larger radii tend to have stronger reduc- 
tion in high velocity particles. This is exempli- 
fied in figure [I] where the upper curves (blue and 
green) show the radial VDF. The open diamonds 
(blue) come from a radial bin in the inner region, 
whereas the stars (green) are from a bin further 
out. The VDF's are normalized such that a com- 
parison is possible, and velocity is normalized to 
the dispersion. This simulated cosmological data 
is from the Local Group simulation of |Moore et al.| 
(2001). The lower curves (red and black) are the 



tangential VDF from the same two bins, however, 
for the tangential VDF there is a striking resem- 
blance, infact to a first approximation these two 
tangential VDF's from different radial bins look 
identical. 

The most frequently used approach to discuss 
DM structures is through the first Jeans equa- 
tion, which relates the velocity dispersions to the 
density profiles. If we were to have some knowl- 
edge about some of the quantities entering the 
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velocity velocity 



Fig. 1. — The velocity distribution function for 2 
different radial bins from a simulated cosmologi- 
cal DM structure ( |Moore et aL|2001| ). The upper 
(green and blue) curves are the radial VDF, and 
the lower (black and red) are the tangential VDF. 
The open diamonds are from the inner bin, and 
stars arc from the outer bin. It is clear that the 
tangential VDF's are very similar to each other, 
whereas the radial VDF's differ in shape both at 
small and large velocities. All figures have velocity 
normalized to the dispersion, and random y-axis 
normalization to enhance visibility. 



Jeans equation, then we can solve for the others. 



One example hereof was presented in Dehnen fc 



McLaughlin (20051, who demonstrated how to de- 
rive a generalized NFW density profile, by both as- 
suming that the pseudo phase-space is a power-law 
in radius ( Taylor fc Navarro|2001 1, and that there 
is a linear relation between the velocity anisotropy 



and the density slope ( Hansen fc Moore] |2006| . 
A somewhat generalized approach was presented 
in Zait et al. (20081, where the authors demon- 
strated how to derive the velocity anisotropy by 
assuming simple forms for both the density profile 
and the pseudo phase-space density. The funda- 
mental problem with this kind of approaches is, 
that any departure from truth in the assumptions 
will lead to departure from correctness in the re- 
sults. Zait et al. ( |2008[ ) demonstrated this in a 
very convincing way, by deriving significantly dif- 
ferent velocity anisotropy profiles by just changing 
between NFW or Sersic density profiles as input. 
Another related problem with these approaches 
is, that the assumption that pseudo phase-space 



Fig. 2. — The velocity distribution function for 2 
different radial bins from the Eddington formula 
for an NFW density profile. The open diamonds 
are from an inner bin (r — 0.1), and stars are 
from an outer bin (r — 10). There is a striking 
resemblance with the radial VDF's from the cos- 
mological simulation in the upper curves in figure 
[l] Same normalization as figure [l] 

is a simple power-law, was recently demonstrated 
to be oversimplified. The unknown question was 
which component (radial, tangential, or something 
else) of the velocity dispersion in the pseudo phase- 
space gives the best approximation to a power-law 



in radius (Hansen et al. 2006 Knollmann et al. 



2008). Schmidt et al. (20081 demonstrated that 



different numerically simulated structures are best 
fitted by different forms of the pseudo phase-space, 
and hence that there is no universal simple behav- 
ior of the pseudo phase-space. 

2.1. The shape of the radial VDF 

For ergodic structures, that is structures where 
the orbits depend only on energy (hence with 
(3 = 0) we can use the Eddington formula to get 



the VDF at any radius ( 


Eddington 


1916), (see 


alsojBinney 


(1982 


); 


Cuddeford 


(1991 


); 


Evans & An 



(2006)). The Eddington formula only depends on 
the radial dependence of the density profile of the 
structure, and by assumption the VDF is the same 
for the radial and tangential directions, f(v,r) = 
function(p(r)). It is natural to interpret this in the 
following way. The structure is in equilibrium, so 
there is detailed balance for each phase-space ele- 
ment. The velocity of each particle is decomposed 
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into the radial and tangential components, and for 
any infinitesimal time step, the radial component 
of any individual particle can tell that it is mov- 
ing in a changing density and changing potential 
(the radial component of any individual particle is 
cither moving directly inwards or outwards). It is 
therefore natural that the radial VDF is imprinted 
by the radial variation in the density and poten- 
tial. For a truncated NFW density profile we can 
e.g. get the VDF at radius 0.1 and 10, in units of 
the characteristic radius, see figure [2j By compar- 
ison of figures [l] and [2] it is clear, that the radial 
VDF of the cosmological simulation indeed looks 
very similar to the VDF from the Eddington for- 
mula. It is therefore tempting to suggest that the 
radial VDF to a first approximation is identical to 
the one which results when applying the Edding- 
ton formula to the given density profile. 

Now, the actual radial VDF (for a given cos- 
mologically simulated structure) will differ slightly 
from the VDF resulting from the Eddington for- 
mula, since the latter was based on the assumption 
that f3 = 0. Specifically, the VDF from the Ed- 
dington formula gives also a r , and to ensure con- 
sistency with the Jeans equation, one must have 
(3 = 0. However, we will here use this VDF as 
a first approximation to the radial VDF, and we 
will present a quantitative comparison in a future 
paper. Recently, Van Hese et al. (20081 showed 



that for a large class of theoretical model, this is 
an excellent approximation (see their figure 5). 

2.2. The shape of the tangential VDF 

It is somewhat less trivial to argue (or claim) 
the shape of the tangential VDF. For an infinites- 
imal time step, the tangential component of any 
individual particle's velocity is moving in constant 
density and constant potential (the tangential ve- 
locity component of any individual particle is mov- 
ing, well, tangentially, and we assume spherical 
symmetry). We still assume that the structure is 
in equilibrium with no time variation. This means, 
that as a first approximation the tangential VDF 
can be thought of as the one resulting from an in- 
finite and homogeneous medium, where both the 
density and potential is constant everywhere. This 
argument is similar to the Jeans swindle, where 
the homogeneous medium implies constant poten- 
tial. Naturally, such an infinite structure is not 
gravitationally stable against perturbations, but 



we can instead approximate it in the following way. 

Let us consider a density profile, which is a 
power-law in radius over many orders of magni- 
tude, and is then truncated. One example hereof 
is an NFW profile, where the central density slope 
is -1, and the corresponding truncation is then 
happening after the scale radius. We are thus 
considering the VDF in a bin at a radius which 
is many orders of magnitude deeper towards the 
center than the scale radius. We can now consider 
a generalized double power-law profile, where the 
central slope can be more shallow than -1, and we 
can use the Eddington formula to extract the VDF 
for any central slope. By lowering the central slope 
towards zero, we get a structure which in principle 
is stable towards perturbations, but at the same 
time is approaching constant density and constant 
potential in the central region. The resulting VDF 
(extrapolated to zero slope) has been discussed by 



Hansen et al. ( 2005 1 and has the shape 



f(v) = n(p) 1 



1-g 
3-<? 



(2) 



with q — 5/3, and n(p) shows that the normaliza- 
tion only depends on the local density. This form 



is known as a q- generalized exponential (Tsallis 



1988). For comparison one should note that a 



simple comparison with polytropes (where f(E) ~ 
^(n-3/2)^ bj-gakg down, since the normal connec- 
tion between density and potential, p ~ ^ r ™, is not 



valid for such shallow slopes (Binncy & Tremaine 



1987) 



When considering the shape of the tangential 
VDF from simulations in figs. [3] and [4] we see 
that this form indeed provides a very good fit for 
all radii, at least for v smaller than roughly 2a. 
The structure in figs. [3] and [4] is from a very non- 
cosmological simulation (the "tangential orbit in- 
stability" of |Hansen et~aT1 ( |2006[ )). The same form 
fits the tangential VDF's from a cosmological sim- 
ulation (lower lines in fig. [T]) equally well. 

Clearly, this form in eq. Q has an extended tail 
of high energy particles, which would not be bound 
by the equilibrated structure. The suppression of 
high energy particles due to the finite radial extend 
of the structure is naturally included through the 
Eddington formula for the radial component, and 
we therefore make the suggestion that the tangen- 
tial VDF must have a high-energy tail which fol- 
lows the radial VDF. Effectively, this means that 
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Fig. 3. — The VDF as function of velocity 
at 3 different radii, upper curves correspond to 
outer radial bin, and lower curves correspond 
to an inner radial bin. The green stars are 
the radial VDF's, whereas the red diamonds are 
the tangential VDF's. The simulation is the 
very non-cosmological "tangential orbit instabil- 



ity" from Hansen et al. ( 2006 ) . The black lines are 
of the form in eq. ([2j). The velocity is normalized 
to the dispersion, and the y-axis has been shifted 
vertically for two of the bins to enhance visibility. 
One clearly sees that the tangential VDF's are vir- 
tually identical, whereas the radial VDF's vary as 
function of radius. 

for large velocities the tangential component of the 
velocity might as well be moving in the radial di- 
rection. This corresponds to the fact that the tan- 
gential velocity component of any individual par- 
ticle actually is moving somewhat radially after a 
finite time-interval. 

When looking at fig. [4] we see that the actual 
suppression is even slightly larger for these high- 
energy particles, however, the difference between 
the suggested and the actual suppression at high 
energy is very small. When looking at the number 
of particles (the integral under the curve in fig. [3]) 
we find, that the difference is virtually zero. 

In conclusion, the tangential VDF is surpris- 
ingly well fit by the phenomenologically predicted 
shape in eq. ([2]), and with a high-energy tail sup- 
pression corresponding to the tail of the radial 
VDF. 

To emphasize the general nature of the shape 
of the tangential VDF, we also present the radial 



Fig. 4. — Same as fig.|3]but with log-scales to make 
the suppression at high-energy more visible. The 
normalization of the y-axis is random, to enhance 
visibility. 



and tangential VDF's of a cosmological simulation 
of a galaxy, including both cooling gas, star- 
formation and stars, as well as supernova feed- 



back from Sommer-Larsen et al. 



(2003); Sommcr- 



Larsen | ( |2006[ ). In fig~[5[ we see that the dark 



matter VDF also in this case has the suggested 
shape. 

3. The velocity anisotropy 

Now, after having established the shape of 
both the radial and tangential VDF's, the velo- 
city anisotropy at any given radius can easily be 
determined, as we will show later in this section, 
since it is just an integral over these distributions, 
cr 2 — J v 2 f(v)dv/ J f(v)dv. The shape of the ra- 
dial VDF changes as a function of radius (section 
2.1), whereas the shape of the tangential VDF is 
virtually constant (section 2.2), and it is therefore 
natural to expect that the velocity anisotropy will 
also change as a function of radius. Since the ra- 
dial VDF generally is more flat-topped than the 
tangential one (see fig. |4j top lines), then we will 
expect (3 to be positive. Only when the density 
slope approaches zero (e.g. a central core) will 
the radial VDF approach the tangential one, and 
hence f3 — > (see fig. [4] bottom lines). 

Let us assume that the radial density profile 
is given, e.g. by a truncated NFW profile. In 
this case the radial VDF is completely determined 
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Fig. 5. — The velocity distribution function as 
function of velocity for a galaxy from a cosmolog- 
ical simulation including both gas and stars from 



Sommcr-Larsen et al. (2003 1; Sommer-Larsen 



(2006). Green diamonds are radial VDF's, red 
stars are tangential VDF's, and the lines are of 
the form in eq. 

through the Eddington formula. The tangential 
VDF is given by eq. ^ and at first sight there are 
2 free parameters, namely the a entering equa- 
tion [2j and then the normalization. For a given a 
we can determine the normalization, since the par- 
ticle number is conserved when integrating over 
the radial or tangential VDF, J f ra ddv = p = 
J ftandv. This leaves us only to determine the 
a entering eq. ([2]). One could argue that it most 
likely is either a r , a tan or 'tot which should en- 
ter here. We will allow ourselves to be guided by 
the results of numerical simulations, and use the 
average a tot, since that gives a fairly good approx- 
imation to all the tangential VDF's from the sim- 
ulations discussed above. We will present a more 
quantitative test of this in a future paper, but the 
effect is modest. E.g. for the truncated NFW pro- 
file at the radius where the slope is —2, we find 
j3 = 0.265 when using a — atou an d if we instead 
use a = oy (atan) we get P — 0.22(0.31). 

It is now straight-forward to calculate the ve- 
locity anisotropy, f3, at any radius for any given 
density profile with no free parameters. Practi- 
cally we do it iteratively in the following way. 1) 
First we find the radial VDF, using the Eddington 
formula. 2) Then we write the tangential VDF, 
which is eq. ^ where we initially use a = a r 



(initially assuming (3 to be very small). 3) Then 
replace this form at high momenta with the radial 
VDF (as described in detail in section 2), and nor- 
malized in such a way that the particle number is 
conserved (between radial and tangential). 4) It 
is now trivial to calculate (3 as an integral over 
these distribution functions, and if this is different 
from the initial assumption, then re-iterate the en- 
tire process with the calculated (3. This means ex- 
plicitely that we return to point 2) with this newly 
calculated (3, and a = o to t- In practice we repeat 
until f3 has converged with accuracy 0.01. 
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Fig. 6. — The velocity anisotropy as function of 
radius. Blue triangles are for an NFW density 
profile, the red diamonds for the density profile 



suggested by Navarro et al. (20041, and the green 
stars is for p(r) ~ 1/(1 + r 2 ) 2 , which has a central 
core (x-axis normalized to the scale radius). The 
squares with error-bars is from the CLEF simula- 



tion ( |Kay et al.||2007| |Springel||2005| ), where the 
67 most relaxed galaxy clusters at z = have 
been selected (x-axis normalized to r 2 5oo)- The 
error-bars correspond to 1 sigma scatter over the 
67 most relaxed clusters (Host et al. 2009). The 



/3-profilcs from pure dark matter simulations (e.g. 
Diemand et al. (2004 1 ) are in good agreement with 
this radial behaviour. 

In fig. [6] we present the radial dependence of (3 
for 3 density profiles, namely an NFW profile with 
a truncation at large radius, a profile like the one 



advocated by Navarro et al. (2004), and finally a 
profile of the form p(r) ~ 1/(1 + r 2 ) 2 , which has a 
central core. We see that the anisotropy increases 
in a way similar to what is observed in numerical 



G 



simulations, namely from something small in the 
central region, to something of the order 0.4 to- 
wards the outer region. The orange squares are 



0.5 FT 



from CLEF numerial simulation ( Kay et aL|2007 
Springel 2005 ) , where the error-bars represent the 



la scatter over the 67 most relaxed clusters (Host 



ct al. 2009). Observations of (3(r) in galaxy clus- 
ters are in excellent agreement with these numeri- 
cal predictions ( Host et al.|2009 ). The radial scale 
of the simulated (3 is r25oo, which does not have 
to coincide with the scale radius of the analytical 
profiles. This gives a free parameter (of the order 
unity) in the normalization of the x-axis, which we 
just put to 1 for simplicity. In a similar way the 
analytical profiles are all normalized to their re- 
spective scale radii, which means that they could 
also have different r25oo- 

Since we have suggested the shape of the full 
VDF's, then we can naturally also get higher order 
moments, such as the kurtosis, as function of ra- 
dius. We thus also predict that the radial profile of 
the higher velocity moments are fully determined 
by the shape of the density profile. 

4. Discussion 

One of the consequences of the above consider- 
ations is, that the appearance of the radial varia- 
tion of (3 is dictated by the density profile. That is, 
given any density profile, the velocity anisotropy 
is entirely determined, as long as the structure 
has had time to equilibrate. We are thus stat- 
ing explicitly, that (3 is unrelated to the infall of 
matter in the outer region, and that the only con- 
nection (3 has to the formation process is through 
the radial structure of the density profile. This is 
naturally supported by the very non-cosmological 
simulations (see figs. [3] [4]) which also produce a 
/3-profile in agreement with cosmological simula- 
tions ( Hansen et al.||200"6| . 

Another consequence is that (3 must always 
be positive in equilibrated structures, since the 
density profile at most can develop a core (see 
fig. |6| . This is also in good agreement with dark 
matter simulations (|Dehnen fc McLaughlin 2005 



Bar nes et al.|2007| |Faltenbacher fc Diemand|2006 



Bellovary et al.| [2008). Virtually no numerical 
simulations find negative velocity anisotropy, and 
when they do, this is usually only in the very in- 
ner region where numerical convergence may be 




Fig. 7. — The velocity anisotropy, (3, as function of 
the density slope, 7. The blue line (triangles) is for 
an NFW profile, and the black solid line (crosses) 
is for a power-law density profile. The orange 
(squares) is for the non-trivial double-bump pro- 
file, the green (stars) is for p = 1/(1 + r 2 ) 2 profile, 
and the red diamonds are for the profile suggested 
Navarro et al. (20041. The dashed straight line 



in 



is the suggestion from Hansen & Stadel (20061. 



These results are roughly fit by (3 = — O.I37. 



questioned. Few analytical treatments have pre- 



dicted a negative (3 in the inner region ( Zait et al. 
2008), however, this result may be an artefact of 



assuming that the pseudo phases-space is a perfect 
power-law in radius, which is generally not correct 
( Schmidt et al. 120081) . If future high-resolution nu- 



merical simulations instead will establish that the 
central velocity anisotropy is negative (in agree- 



ment with the predictions of Zait et al. (2008)), 



then that would be a proof that the present anal- 
ysis is flawed somehow. 

It has previously been suggested that a con- 
nection between the anisotropy and the slope of 
the density profile should exist. This connec- 
tion appears to hold even for structures which 
have profiles with non-trivial radial variation in 
dlogp/dlogr ( |Hansen fc Moore|2 006). We can now 
test this connection. In figure [7] we show f3 as 
function of the density slope for the NFW profile 
(solid line, blue triangles), the profile suggested 
by Navarro et al. (2004) (red diamonds), and also 



for a "double-bump" structure (the sum of two 
spatially separated profiles of the form 1/(1 + r) 3 , 
orange squares)). This double- hump profile has 
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a very non-trivial radial variation of dlnp/dlnr, 
which cannot be well approximated by any gen- 
eralized double power-law profile. All these struc- 
tures appear to land near a connection roughly 
given by (3 = —0.1 37. We also show the results 
for the p — 1/(1 + r 2 ) 2 profile (green stars), as 
well as for single power-law profiles (fat black line, 
crosses). The dashes line is = —0.2(7 — 0.8) as 
suggested in Hansen & Stadcl (2006), based on a 
set of cosmological and non-cosmological simula- 
tions. These two results differ by approximately 
0.1 in (3. We indeed see that all structures land 
in a relatively narrow band in the 7-/? plane, and 
hence likely explaining the origin of the 7-/? rela- 
tions. 

The most important practical implication of 
this suggestion is, that it will allow us to close the 
Jeans equation. As is well-known, the Jeans equa- 
tion depends on the density, dispersion, anisotropy 
and the total mass. Now, having demonstrated (or 
at least suggested strongly) that the anisotropy is 
uniquely determined once the density is known, we 
see that it is possible to close the Jeans equation 
for systems that are fully relaxed. 

We have been making simplifying assumptions 
above, which all need to be tested through high 
resolution simulations. First, we assume that the 
radial VDF is very similar to the one appearing 
from the Eddington formula, even in the presence 
of a non-zero (3. We also assume that the a en- 
tering eq. ([2]) is the total one. If the correct a to 
use is instead closer to a ta m then (3 will be slightly 
larger, but the radial variation will remain. From 
these assumptions we estimate the accuracy of the 
present work to be about 0.1 or up to about 30% 
in (3{r). 

One could naturally ask why and how the ra- 
dial and tangential VDF's get their shapes? It is 
slightly disappointing that it is not a deep phys- 
ical principle, like a generalized entropy, which 
is responsible. Instead it is simply the density 
profile (cither the radially varying, or the tan- 
gcntially constant) which through the Eddington 
formula demands that the VDF's take on these 
forms. These forms will therefore appear when 
there is sufficient amount of violent relaxation to 
allow enough energy exchange between the parti- 
cles. 



5. Summary 

The velocity of any particle can be decomposed 
into the radial and tangential components, and 
when summing over all particles in a radial bin, 
we get the particle velocity distribution function, 
the VDF. We suggest that both the radial and 
tangential VDF's are given through the Eddington 
formula. The radial one comes from the radially 
changing density profile, and the tangential VDF 
arises when considering a structure with constant 
density and potential. This is because the tangen- 
tial component of the velocity as a first approx- 
imation is moving in constant density and con- 
stant potential. In addition the tangential VDF 
is reduced for high-energy particles in accordance 
with the radial VDF, to ensure that the particles 
remain bound to the structure. These phenomeno- 
logical predictions are in remarkably good agree- 
ment with the results from numerical simulations 
of collisionless particles, both of structures of cos- 
mological origin as well as highly non-cosmological 
origin. 

Under these suggestions it is straight forward 
to derive the velocity anisotropy profile, (3(r), with 
no free parameters. This is shown to increase ra- 
dially from something small (possibly zero) in the 
center, to something large and positive (possibly 
around 0.4) towards the outer region. 

We have thus demonstrated that the velocity 
anisotropy is entirely determined from the density 
profile. This allows us to close the Jeans equation, 
since (3 is no-longer a free parameter. 
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